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AN  EMPIRICAL  PREDICTION  ALGORITHM  FOR  LOW-FREQUENCY 
ACOUSTIC  SURFACE  SCATTERING  STRENGTHS 


1.  INTRODUCTION 

Acoustic  models  that  attempt  to  predict  reverberation  levels  for  a  particular  set  of  environ¬ 
mental  and  acoustic  parameters  require  a  model  for  ihe  level  of  surface  backscatter.  For  most 
low-frequency  (i.e.,  below  1  kHz)  models,  as  well  as  for  general  sonar  equation  calculations,  the 
surface  backscatter  component  of  the  general  reverberation  field  is  modeled  using  the  empirically- 
derived  Chapman-Harris  formula  (Chapman  and  Harris  1962).  This  formula  is 

55  =  3.3/3  log( 0/30)  -  42.4  log/3  +  2.6, 

where 

and  55  is  the  surface  scattering  strength  in  dB,*  6  is  the  grazing  angle  in  degrees,  U  is  the  wind 
speed  in  knots,  and  /  is  the  acoustic  frequency  in  Hz.  This  relationship  was  derived  from  data 
collected  in  deep  water  off  the  coast  of  Bermuda  by  using  signal,  underwater  sound  (SUS)  explosive 
charges  as  sources  and  a  single  receiving  hydrophone.  The  data  were  collected  during  a  single  52-h 
period  when  surface  conditions  ranged  from  sea  state  0  to  sea  state  6,  and  wind  speeds  ranged 
from  0  to  30  knots.  The  data  wore  measured  in  octave  frequency  bands  from  400  to  6400  Hz. 
Later  work  by  Chapman  and  Scott  (1964)  extended  the  data  range  to  lower  frequencies. 

While  the  Chapman-Harris  empirical  formula  adequately  describes  the  levels  of  surface  scat¬ 
tering  for  some  combinations  of  frequency  and  surface  conditions,  recent  work  in  programs  such 
as  the  Critical  Sea  Test^  (CST)  shows  that  there  are  other  frequency/environment  regimes  in 
which  the  Chapman-Harris  predictions  are  not  adequate.  Here  we  propose  a  generalization  of  the 
Chapman-Harris  formula  for  modeling  surface  scattering  strength  as  a  function  only  of  frequency, 
wind  speed,  and  grazing  angle,  ba.se<l  on  new  observations  of  surface  scattering  made  during  five 
CST  at-sea  exercises.  The  data  on  which  this  generalization  are  based  will  not  be  presented  here, 
as  they  have  been  reported  elsewhere  (Ogden  1992;  Erskine  et  al.  1992).  Hather,  this  report 
summarizes  the  results  obtained  in  the  CS'l  exercises  and  discusses  the  algorithm  based  on  those 
results. 

This  algorithm  is  intended  to  be  an  initial  step  between  the  existing  Chapman-Harris  formula 
and  a  more  detailed  description  of  surface  scattering  that  is  the  subject  of  continuing  research.  The 
algorithm  provides  a  better  estimateof  .scattering  strength  in  the  .60  to  1000  Hz  band  for  the  grazing 
angles  and  wind  speeds  covered  in  this  rejiort  than  does  Chaimian-Harris.  However,  additional 


Manivscript  approved  February  23,  1992. 

' SS  is  a  dimensionless  quantity,  usin^  ihe  definition  of  iiV  =  10  log  gjven  bv  Frirk  (1983). 
'Sponsored  by  Spare  and  Naval  Warfare  Systems  Command,  PMW-183. 
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work  is  needed  on  the  surface  scattering  problem  to  arrive  at  a  more  detailed  umb'rstaiidiiig  of  the 
dependence  of  surface  reverberation  on  environment;d  and  acoustic  parameters. 

There  are  several  qualifications  that  must  be  applied  to  the  restill  propo.sed  in  this  report. 
The  most  important  qualification  is  that  it  ha.s  been  recognized  from  (’S  f  and  other  work  on  surface 
scattering  that  a  prediction  of  scattering  strength  based  on  a  single  environnu'ntal  parannUer  (i.<'.. 
on  instantaneous  wind  speed  alone)  is  not  aderpiate.  For  exaniph',  it  has  bc'en  foutul  that  a 
[tarameter  like  sea  state  that  contains  some  wind  history  appears  to  be  a  lu'ttc'r  predictor  of 
surface  scattering  strength  than  wind  speed  (Ogden  199‘2:  I'rskine  et  al.  1992).  For  this  reason,  it 
would  be  preferable  to  use  sea  state  as  a  parameter  rather  than  wind  spc'c'd  in  modcds  of  surface 
reverberation.  However,  to  make  the  algorithm  proj)osed  in  this  report  as  compatible'  tes  jiossiblc' 
with  existing  models,  the  deimndence  of  scattc'ring  stretigth  on  wind  sju'c'd  alone  has  bc'c'n  rctaiuecl. 
\  .second  qualification  is  that  there  has  been  no  attempt  made  in  the  dc'vc'loptnent  of  the-  proposc'd 
algorithm  to  do  any  detailed  functional  fitting  of  the  (.’SI'  surface  scattc'r  data,  fhe  approac  h 
we  have  taken  is  to  identify  empirical  regime's  in  the  range  of  frecjuc'ncy/surfacc'  conditions  wiu'rc' 
e.xisting  theoretical  or  empirical  formulas  seem  to  be  appropriate  and  to  merge'  thc'in  smoothly  into 
a  single  predictive  algorithm. 


2.  DESCRIPTION  OF  CST  DATA 


One  experiment  that  has  been  a  component  of  each  CS'f  at-sea  exc'rcisc'  has  bc-en  tin'  me;i- 
surement  of  surface  backscattering  strengths  using  SUS  charge's  as  source's  atid  a  towc'd  horizontal 
line  array  as  a  receiver.  In  this  experiment,  the  SFS  charge's  are'  dc'totiatc'd  dirc'ctly  bc'uc'ath  the' 
receiver,  and  the  backscattered  surface  returns  are  re'ce'i',ed  on  a  set  of  hydrophone's.  Thc'  signals 
from  these  hydrophetnes  are  thc'n  beamforinc'd  and  Fouric'r  analyzc'd  into  narrowband  (gc'uc'r.dly  ! 
Hz)  reverberation  levc'ls.  Scattering  strc'iiglhs  arc'  the'ii  calc’idalc'd  from  tlu'sc'  revc'rbc'rat ion  Icvc'ls 
by  combining  them  with  the  source  level  for  the'  SFS  charge',  the'  cc)m[)utecl  t ratismissiou  ho"  irom 
the  se^iirce  to  the  surface  and  from  the'  surface'  to  the'  rc'e'e'ive'r.  ami  tin'  area  insonilied  by  t  ,  -  signal 
during  a  particular  period  of  time.  Hc'cause'  SFSs  arc'  broadband  souicc's.  scattc'iinz  stii'iigth' 
may  be  measured  simultanc'ously  for  a  wide'  rang*'  of  frc'ejuc'ncic's.  In  tin'  ( 'S  1  prog  .  .n.  tin'  SI  S 
measure'ments  are'  gc'uc'rally  analyzed  ovc-r  a  range'  of  frecpic'iicic's  from  '-7(1  Hz  '  ■  about  1  kHz. 
riie  grazing  angle's  covc'rc'cl  bv  this  e'xpc'rinn'iit  are'  usually  be'tvve'e'ii  about  20'  an  i  arouml  b  to  7  . 
riie  lower  limit  is  ge'in'rally  ch't  ermined  by  tin'  time'  of  I  he'  arrival  ol  tin'  bot  tom  la  t  hoiuet  oi  ot  u  rn . 
A  detailed  eh'scriiction  of  the  acejuisilion  anci  analysis  jtroce'ss  may  be-  louud  olscw  la'ie  ((tgdi  u  ami 
Krskine'  19.S9). 


During  the'  first  four  t'S  l  at-s<'a  exe'rcisc's,  a  total  of  2 1  SI  ''  -urfacc-  scatti'r  cbita  -'  t- 
collecte'd.  Kach  data  set  took  roughly  an  hour  ancl  involvc'd  ■ne  dc'jihiyunut  of  10  to  pit  SI  !s 
charges.  I  hese'  rne'asu reme'iit s  wc're'  carrie'd  out  uiicIct  a  range'  c  f  envi roti tm-iO al  c  cunlit  loio;  during 
the  various  SFS  Ic'sts.  wind  sjcc'c'cls  varie'd  from  b  to  2S  knots  ( c  c irrei  t I'd  tn  a  In'i'ght  ot  I't  .i  ni’ i. 
with  sea  states  from  1  to  t,.'), 

*An  oftrn-ii.snl  referetu  c  luiKlit  tlcat  conn's  from  tlic  w  vc  s|>(-f|rii  li.rta  <oi|iili2  la  i  .o"! 

analyzed  by  Pierson  and  Moskowitz  (ittSI) 
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Ea/ch  of  the  individual  runs  was  analyzed  to  give  scattering  strength  as  a  function  of  grazing 
angle,  frequency,  and  environmental  conditions.  The  results  of  these  runs  were  then  compared 
to  empirical  and  theoretical  predictions  of  scattering  strengths  for  the  conditions  that  occurred 
during  the  runs.  Figure  1  gives  the  results  of  these  comparisons  and  summarizes  many  of  the  SUS 
runs  by  dividing  the  appropriate  frequency  and  wind  speed  space  into  three  regimes.  (It  is  worth 
noting  again  that  the  use  of  wind  speed  as  an  environmental  descriptor  here  is  not  an  optimum 
choice;  it  does  not  include  the  other  effects,  such  as  wind  history,  that  are  known  to  be  important 
in  the  prediction  of  scattering  strengths.) 

In  the  regime  that  encompasses  the  lowest  wind  speeds  at  all  frequencies  and  increasingly 
higher  wind  speeds  at  lower  frequencies,  the  scattering  strength  results  were  found  to  be  in  rea¬ 
sonable  agreement  with  the  predictions  of  air-sea  interface  scattering  as  described  by  perturbation 
theory  (see,  for  example,  Thorsos  1990).  The  principal  predictions  of  perturbation  theory  for  the 
range  of  frequencies,  grazing  angles,  and  wind  speeds  relevant  here  are  that  the  scattering  strengths 
should  have  a  tan'*^  dependence  on  grazing  angle  and  relatively  bttle  dependence  on  frequency 
and  wind  speed. 

At  higher  wind  speeds  and  higher  frequencies,  the  SUS  results  were  found  to  be  in  reasonable 
agreement  with  the  Chapman-Harris  empirical  formula  given  previously.  While  the  grazing  angle 
dependence  of  Chapman-Harris  is  roughly  similar  to  that  of  perturbation  theory,  the  Chapman- 
Harris  formula  has  considerably  more  wind  speed  and  frequency  dependence  than  perturbation 
theory  predictions.  This  enhanced  dependence  results  in  scattering  strength  levels  that  are  con¬ 
siderably  higher  than  can  be  expected  from  scattering  from  the  air-water  interface,  which  strongly 
suggests  that  the  principal  scattering  mechanism  must  be  different  from  just  the  scattering  off  a 
rough  surface.  It  is  widely  believed  that  subsurface  bubble  clouds  and/or  plumes  give  rise  to  this 
enhanced  acoustic  scattering. 

In  between  these  two  regions  is  a  transitional  region  where  the  two  mechanisms  are  (pre¬ 
sumably)  competing  as  the  dominant  source  of  scattering.  The  scattering  strengths  in  this  regime 
are  generally  somewhere  between  the  predictions  of  perturbation  theory  and  Chapman-Harris,  but 
the  exact  levels  are  found  to  depend  on  environmental  ajid  acoustic  parameters  in  a  manner  that 
is  not  yet  well  understood. 

As  previously  mentioned,  the  regimes  shown  in  Fig.  1  were  determined  by  examining  many 
individual  scattering  strength  vs  grazing  angle  curves  and  making  a  judgment  as  to  whether 
they  were  best  described  by  Chapman-Harris  or  perturbation  theory,  or  were  in-between  these 
descriptions.  The  boundaries  of  the  three  regimes  were  deterrmned  through  the  production  of  the 
plot  shown  in  Fig.  2.  In  this  plot,  each  point  represents  a  grazing  angle  vs  scattering  strength  curve 
at  a  single  frequency,  with  the  letter  indicating  the  best  description  of  the  entire  curve.  (Note  that 
in  some  parts  of  the  frequency/wind  speed  regime,  Chapman-Harris  and  perturbation  theory  give 
similar  predictions.  These  points  were  normally  assigned  to  perturbation  theory.)  After  plotting 
all  the  points,  the  boundaries  of  the  regimes  were  determined  visually  rather  than  by  using  an 
analytical  scheme. 


3 


Wind  Speed  (knots) 


OGDEN  AND  ERSKINE 


No  CST  SUS  data  available  ^x>ve  28  knots 


ONSET  OF 
SIGNIFICANT 
WHITECAP 
CONCENTRATION 

(TRANSITION 

FROM 

SEA  STATE  2  TO 
SEA  STATE  3) 


100  200  300  400  500  600  700  800  900  1000 


Frequency  (Hz) 


NO  WIND  SPEED 
DEPENDENCE 

VERY  SLIGHT 
FREQUENCY 
DEPENDENCE 

GENERALLY 
CONSISTENT  WITH 
OR  SLIGHTLY  ABOVE 
PERTURBATION 
THEORY 


DEPENDENCE  ON 
WIND  SPEED  AND 
FREQUENCY  NOT 
WELL  UNDERSTOOD 

TRANSITION  REGION 
FROM  AIR-WATER 
INTERFACE 
SCATTERING  TO 
ANOTHER  MECHANISM 


•WIND  SPEED  AND 
FREQUENCY 
DEPENDENCE 
SIMILAR  TO 
CHAPMAN-HARRIS 

■ANOTHER 
SCATTERING 
MECHANISM 
REQUIRED  (POSSIBLY 
BUBBLE  PLUMES) 


Surface  scatter  results  from  the  CST  SUS  tests 


4 


NRL  REPORT  9377 


No  CST  SUS  data  available  above  28  knots 


Frequency  (Hz). 
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perturbation  theory 
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empirical  formula 


Fig.  2  -  Individual  points  used  to  produce  Fig.  1 
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3.  SCATTERING  STRENGTH  ALGORITHM 

The  surface  scattering  algorithm  presented  in  this  report  is  a  two-step  procedure.  First,  the 
wind  speed  and  frequency  are  used  to  decide  which  of  the  regimes  shown  in  Fig.  1  is  appropriate. 
Second,  the  scattering  strength  is  calculated  from  the  wind  speed,  frequency,  and  grazing  angle 
by  using  one  of  three  methods  that  correspond  to  the  three  regimes.  The  three  calculations  of 
scattering  strength  used  in  this  algorithm  are  (1)  first-order  perturbation  theory,  (2)  the  Chapman- 
Harris  empirical  formulation,  and  (3)  an  interpolation  between  (1)  and  (2). 

The  first  step  in  translating  the  regime  plot  in  Fig.  1  into  an  algorithm  was  to  put  the  bound¬ 
aries  of  the  regions  into  analytical  form.  The  lower  boundary,  between  the  perturbation  theory 
region  and  the  transition  region,  was  approximated  by  two  line  segments.  The  upper  boundary, 
between  the  transition  region  and  the  Chapman-Harris  region,  was  fit  by  a  cubic  polyinonial. 
These  two  boundary  curves  intersected  at  about  .35  knots  and  50  Hz,  so  50  Hz  was  taken  as  the 
low-frequency  cutoff  to  the  algorithm.  We  have  analyzed  data  up  to  ~1  kHz,  so  this  was  used  as 
the  upper  frequency  cutoff.  As  noted  previously,  the  range  of  grazing  angles  over  which  we  have 
data  is  about  5°  to  30°,  but  we  have  extended  the  upper  and  lower  limits  of  the  algorithm  to 
1°  and  40°  respectively,  as  these  limited  extrapolations  are  probably  valid.  (Perturbation  tln'ory 
becomes  a  poor  approximation  at  high  grazing  angles  unless  other  effects  are  added,  thus  we  have 
used  a  40°  cutoff.)  An  arbitrary  upper  limit  of  40  knots  was  included,  while  a  lower  wind  speed 
limit  of  5  knots  was  included  for  reasons  that  will  be  discussed  later. 


For  lower-wind,  lower-frequency  cases,  the  algorithm  uses  two-dimensional  first-order  per¬ 
turbation  theory  for  calculating  scattering  strengths.  A  presentation  of  the  relevant  equations  is 
given  by  Thorsos  (1990).  Using  the  notation  of  Thorsos,  the  2-D  backscattering  cross  section 
is 


where  k,^  is  the  incident  acoustic  wavenumber  in  the  z  direction.  Vf'(A’)  is  the  2-1)  roughness 
spectral  density,  K  is  the  “transverse”  surface  wave  vector  (A'  =  Kj-t  -f  Kyi)),  and  the  scattering 
strength  is  given  by  10  log  crjp,  where  the  superscript  on  rr  refers  to  t  he  nature  of  the  terms  included 
in  a  first-order  perturbation  calculation.  If  w'e  make  the  assumption  that  th('  wind  direction 
along  the  x  axis,  and  we  use  the  definition  of  H''(  A')  given  by  Thor.sos,  the  expression  for  a  becomes 


,(2) 

^2D 


where  $  is  the  azimuthal  dependence  of  the  wave  spectrum,  and  we  have  chosen  to  take  S('2k,j-) 
as  the  Pierson- Moskowitz  1-D  wave  spectral  density  (Pierson  and  Moskowitz  1964)  given  by 


S{2k,,) 


a 

32]^ 


exp 


where  a  =  8.10  x  10~^,  fi  =  0.74,  g  -  9.81  m/s^,  and  U  is  the  wind  speed  in  m/s  at  a  height 
of  19.5  m.  Since  the  input  to  the  algorithm  is  the  acoustic  frequency  /,  we  use  the  definitions 
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sin  9  and  fc,x  =  ^  cos  6  in  the  above  expression  to  get 

1.01  X  10^  1 
pU*  cos^  0 

For  the  surface  scatter  algorithm  we  have  assumed  that  surface  roughness  is  isotropic,  in  which 
case  ^  and  we  finally  arrive  at 


(T20  =  1.01  X  10  ^tan^dexp 


a^2D  =  1.61  X  10  '‘tan'*^exp 

As  noted  previously,  10  log  <7^^  is  the  perturbation  theory  expression  contained  in  the  surface 
scatter  algorithm. 


.01  X 

/2l74^os2  9 


The  only  additional  modification  to  the  actual  calculation  is  that  the  wind  speed  lower  limit 
is  set  to  5  knots,  with  values  lower  than  this  rounded  to  5  knots.  Since  the  Pierson-Moskowitz 
formula  was  determined  by  measuring  surface  wave  spectra  for  wind  speeds  between  20  and  40 
knots,  it  is  not  at  all  clear  that  it  is  applicable  to  very  low  wind  speeds.  Since  there  is  evidence 
that  surface  scattering  strengths  are  still  reasonably  web-characterized  by  perturbation  theory  for 
wind  speeds  around  4  to  5  knots  (Ogden  and  Erskine,  to  be  published),  we  have  chosen  to  put 
a  lower  limit  on  the  wind  speed  to  avoid  the  return  of  unreasonable  values  at  low  frequencies. 
As  with  the  grazing  angle  cutoffs,  the  low  wind  speed  cutoff  is  unlikely  to  affect  any  serious 
reverberation  calculations,  since  surface  scatter  from  wind  speeds  below  5  knots  is  unlikely  to  be 
of  much  importance  compared  to  other  noise  sources. 


The  remainder  of  the  algorithm  is  straightforward  to  describe.  At  higher  frequencies  and 
higher  wind  speeds,  the  Chapman-Harris  empirical  formtda  presented  earlier  is  used.  In  the  tran¬ 
sition  region,  both  the  perturbation  theory  scattering  strength  SSp^rt  and  the  Chapman-Harris 
scattering  strength  SScu  are  calculated.  Then  a  simple  interpolation  on  the  resulting  dB  values 
is  performed: 

SS total  -  OtSScH  +  (1  -  0')5’5'per(, 

where 

input  wind  speed  —  wind  speed  at  pert,  theory  boundary 
•nput  speed  at  CH  boundary  —  wind  speed  at  pert,  theory  boundary’ 

where  the  boundary  wind  speeds  are  evaluated  at  the  input  frequency.  Thus  at  a  particular 
frequency,  a  wind  speed  close  to  the  perturbation  theory  regime  boundary  will  give  a  scattering 
strength  close  to  the  perturbation  theory  level,  while  a  wind  speed  just  below  the  Chapman-Harris 
boundary  will  result  in  a  scattering  strength  that  is  close  to  Chapman-Harris.  While  this  approach 
is  clearly  an  approximation,  it  does  reflect  the  observed  fact  that  for  conditions  in  the  transition 
region,  the  scattering  strength  values  are  usually  somewhere  between  the  two  extremes. 


To  test  the  algorithm,  we  generated  a  series  of  plots  that  hold  one  of  the  three  input  param¬ 
eters  (frequency,  grazing  angle,  wind  speed)  constant  and  show  how  scattering  strength  varies  as  a 
function  of  the  other  two.  Examples  of  these  plots  are  shown  in  Figs.  3  through  8.  In  Figs.  3  and 
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Fig.  8 


Surface  scattering  strengths  as  a  function  of  grazing  angle  and  frequency  for  a  wind  speed  of  5  knots 
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Surface  scattering  strengths  as  a  function  of  grazing  angle  and  frequency  for  a  wind  speed  of  30  knots 
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4,  the  grazing  angle  is  held  constant  at  10°  and  30°,  respectively,  while  wind  speed  and  frequency 
vary.  In  Figs.  5  and  6,  the  frequency  is  held  constant  at  100  and  900  Hz,  respectively,  while  in 
Figs.  7  and  8,  the  wind  speed  is  held  constant  at  5  and  30  knots,  respectively. 

Comments  about  the  behavior  of  the  algorithm  at  some  of  the  “edges”  of  the  parameter 
space  are  in  order.  The  effect  of  the  5-knot  wind  speed  cutoff  can  be  seen  clearly  in  Figs.  3  and 
4,  while  the  1°  grazing  angle  cutoff  is  apparent  in  Figs.  5,  6,  and  8.  In  Fig.  3.  the  transitions 
between  perturbation  theory  and  Chapman-Harris  may  seem  too  abrupt,  but  they  in  fact  reflect 
the  behavior  of  the  data.  The  exceptions  to  this  are  the  transitions  at  the  lowest  frequencies  (not 
shown  on  Fig.  3)  that  occur  above  30  knots.  These  transitions  probably  are  too  abrupt,  but 
unfortunately  there  are  no  data  on  which  to  base  a  better  prediction.  The  rolloff  of  the  50  Hz 
curves  in  Figs.  3  and  4  at  low  wind  speeds  may  not  be  realistic  in  the  presence  of  significant 
swell  components,  but  again  there  are  no  data  to  support  a  better  model.  In  any  case,  the  rolloff 
probably  gives  a  generally  correct  feeling  for  the  importance  of  surface  scatter  in  these  conditions. 
In  Fig.  7,  the  .50  Hz  curve  is  off  the  bottom  of  the  plot,  and  again  this  may  be  too  low.  The  data  on 
w’hich  the  algorithm  is  based  were  analyzed  with  70  Hz  as  the  lowest  frequency.  While  the  70  Hz 
levels  from  the  algorithm  generally  give  results  in  agreement  with  the  data,  the  presence  of  swells 
may  in  fact  keep  the  50  Hz  scattering  levels  in  the  deep  ocean  higher  than  a  Pierson-Moskowitz 
surface  spectrum  would  suggest.  Again,  there  are  no  data  on  which  to  base  a  better  prediction. 

The  algorithm  itself  has  been  written  in  two  forms:  as  a  FORTRAN  function  and  as  a 
Microsoft  Excel®macro  function.  Figure  9  shows  the  FORTRAN  code  listings  and  Fig.  10  shows 
the  Excel  macro. 


4.  SUMMARY 

An  algorithm  has  been  proposed  for  calculating  surface  backscattering  strength  for  grazing 
angles  between  0°  and  40°.  for  wind  speeds  between  0  and  40  knots,  and  for  frequencies  between 
50  and  1000  Hz.  flie  algorithm  is  intended  as  a  replacement  for  the  Chapman-Harris  empirical 
formula  presently  in  use  in  virtually  all  low-fre(juency  acoustic  r<’V('rb('ration  models. 

fhe  algorithm  is  based  on  data  collerte<l  during  the  C'ST  series  of  at-sea  experiments.  In  ('ach 
f)f  the  first  five  (\S'i  exercises,  experinnuits  were  conducted  to  nu'asure  surface  scattering  strengths 
using  broafiband  SI'S  charges  as  sources,  'fhe  results  of  these  exi)«’riments  were  compared  to 
two  existing  methods  for  predicting  surface  scattering  strengths:  air-water  interface  scattering 
theory  (in  the  form  of  first-order  perturbation  theory),  and  the  Chapman-Harris  emirirical  formula. 
Environmental  conditions  where  either  of  the  two  predictive  methods  appear  to  be  appropriate 
were  identified,  as  well  as  a  region  that  had  levels  that  were  intermediate  between  the  two.  I'he 
algorithm  uses  the  CS  T  data  to  identify  where  the  boundaries  of  the  scattering  regions  are  and 
provides  a  smooth  transition  between  the  dilferent  calculation  nu'thods. 

The  algorithm  described  here  should  be  regarded  as  an  interim  solution  to  the  problem  of 
predicting  surface  scattering  strengths.  On  the  one  hand,  the  evidence  is  clear  that  the  Chapman- 
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RSAL* *4  FUNCTION  SURFACB_SS 

>  (windspeed , f  requency , grazing_angle) 

« 

*  This  function  calculates  surface  scattering  strengths,  given 

*  wind  speed  in  knots, 

*  frequency  in  Hz, 

*  grazing  angle  in  degrees. 

* 

•  If  the  wind  speed  <  5.0  knots,  the  surface  scattering  strength 

•  will  be  computed  for  a  wind  speed  of  5.0  knots. 

•  If  the  wind  speed  >  40.0  knots,  the  surface  scattering  strength 

•  will  be  computed  for  a  wind  speed  of  40.0  knots. 

• 

»  If  the  frequency  is  <  50.0  Hz  or  >  1000  Hz,  an  ERROR  CONDITION 

•  is  returned.  The  scattering  strength  will  be  set  to  +1000. 

♦ 

*  If  the  grazing  angle  is  <  1.0,  the  surface  scattering  strength 

*  will  be  computed  for  a  grazing  angle  of  1.0  degrees. 

‘  If  the  grazing  angle  is  >  40.0,  the  surface  scattering  strength 

*  will  be  computed  for  a  grazing  angle  of  40.0  degrees. 


real*4 

real*4 

real*4 

real*4 

real*4 

real*4 

real*4 

real*4 

real*4 


beta 

chapman_harris_liniit 

frequency 

grazing_angle , grazing_angle_r 
interpoIation_f  actor 
perturbation_limit 
PI 

tempi , tcmp2 , temp3 , temp^ss 
wind , wind_speed 


parameter  (PI  =  3.1415927) 


• —  if  wind  speed  is  out  of  range,  reset  it  to  the  extremes: 
if  (wind_speed  .It.  5.0)  then 
wind  =  5.0 

else  if  (wind_specd  .gt.  40.0)  then 
wind  =  40.0 
else 

wind  =  wind_speed 
endif 


• -  if  frequency  is  out  of  range,  return  with  an  ERROR  condition: 

if  ((frequency  .It.  50.0)  .or.  (frequency  gt.  1000.0))  then 
surface_ss  =  1000.0 
return 
endif 


‘ —  if  the  grazing  angle  is  out  of  range,  reset  it  to  the  extremes 
•  —  also;  convert  to  radians: 

Fig.  9  -  FORTRAN  listing  of  surface  scatter  algorithm 
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if  (gruing_aacl«  .It.  1.0)  then 
gruing_»ngl«_r  =  PI/180.0 
else  if  (gruiiig_aiigle  .gt.  40.0)  then 
grazing_angle_r  =  40.0  •  PI/180, 
else 

grazing_angle_r  =  grazing_angle  *  PI/180, 
endif 


determine  limit  of  perturbation  theory  range: 
if  (frequency  .ge.  240.0)  then 
perturbation_limit  =  14.0 
else 

perturbation_limit  =  41.71  -  0. 1 155»f requency 
endif 


check  if  perturbation  theory  applies; 
if  (wind  .le.  perturbation_limit)  then 
use  perturbation  theory: 

tempi  =  1 . 43E7/ (frequency»*2.0) / (wind»»4 . 0) 

>  /  (cos(grazing_angle_r) **2) 

isotropic  form  of  perturbation  theory; 
temp2  =  1.61B-4  *  tan(grazing_angle_r) **4. 

>  ♦  exp (-min (tempi ,200.0) ) 

perturbation  theory  scattering  strength: 
if  (temp2  .gt.  l.OE-20)  then 

surface_ss  =  10. *logl0(temp2) 
else 

surface_ss  =  -200.0 
endif 

else 

find  the  chapman  harris  limit  at  this  frequency: 
chapman_harris_limit  =  39.07-0.066»frequency 

>  *  7.06E-5  •  (frequency»»2. )  -  2.58E-8  *  (frequency»*3. ) 

apply  chapman  harris  formula; 

beta  =  158 .* (wind* (frequency**0 , 3333) )*» (-0. 58) 
tempi  =  beta  *3.3  *  logl0(grazing_angle/30. ) 
temp2  =  42.4*logl0(beta) 
chapman  harris  formula; 
temp3  =  tempi  -  temp2  *  2.6 

check  if  chapman  harris  formula  applies  directly: 
if  (wind  .ge.  chapman_harris_limit)  then 
surfacess  =  max (temp3, -200. 0) 

else 

find  interpolation  factor; 

interpolation_f actor  =  (wind  -  perturbation_limit) 

>  /  (chapfflan_harris_limit-perturbation_limit) 

tempi  =  1.43E7  /  (f requency**2 . ) / (wind**4 . ) 

>  /  (cos(grazing_angle_r)**2.) 

isotropic  form  of  perturbation  theory: 
temp2  =  1.61E-4  *  tan(grazing_angie_r) **4 . 

>  *  exp(-min(templ, 200.0)) 

perturbation  theory  scattering  strength: 
teMp_ss  =  10. •Iogl0(t««p2) 
interpolation  done  in  dB  space: 

Fig.  9  -  FORTRAN  listing  of  surface  scatter  algorithm,  continued 
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=  te*p_ss  ♦  interpolatioii_f actor  *  (tempS-teap  ss) 
surfacess  =  iiax(teap_ss, -200.0) 
endif 
endif 

return 

end 


Fig.  9  -  FORTRAN  listing  of  surf«ice  scatter  algorithm,  continued 


Surface_SS 

Surface  scatter  algorithm  based  on  CST  SUS 

»ARGUMENTCwind*) 

wind  speed  in  kts 

=ARGUMENTCfreq-) 

frequency  in  Hz 

=ARGUMENT(“anqle“) 

grazing  angle  in  degrees 

«IF(OR(freq<50,freq>1 000), RETURN!  1000)) 

test  for  in-range  values 

-IF(angle<=1,1,anqle) 

lower  angle  limit  of  1  deg 

-IF(A60>40,40,A60) 

upper  angle  limit  of  40  deg 

-IF(wind>40,40,wind) 

upper  wind  limit  of  40  kts 

-1F(A62<=5.5,A62) 

lower  wind  limit  of  5  kts 

-IF(ffeq>=240,14.4l.71-0.1 155‘ffeq) 

idetermine  limit  of  pert,  theory  range 

«1F(A63<=A64,G0T0(A78)) 

Ibranch  if  pert,  theory  applies 

-39.07  0.066*freq+0.0000706Mreq'2  0.0000000258*freq"3 

CH  limit  at  this  frequency 

=IF{A63>=A66,GOTO(A82)) 

branch  it  CH  applies 

=  (A63A64)/(A66-A64) 

interpolation  factor 

-158*(A63*freq"(1/3))'(-0.58) 

beta 

=.3.3'A69’LOG10(A61/30) 

=42.4*LOG10(A69) 

=  A70^71  +2.6 

Chapman-Hams  formula 

=  14300000/freq''2/A63"4/COS(A6rPI()/i80)'2 

exponential  of  PM  -  note  conversion  to  kts 

^.0001 61  *TAN(A6rPIQ/180)''4-EXP(  MINI A73. 200)) 

isotropic  form  of  pert  theory 

=  10*LOG10(A74) 

pert,  theory  SS 

=  A75  +  A68'(A72-A75) 

interpolation  done  in  dB  space 

=  RETURN(MAX(A76,-200)} 

return  interpolated  SS 

=  14300000/freq''2/A63''4/COS(A6I  •Pin.'180)'2 

exponential  of  PM  -  note  conversion  lo  kts 

=  0.000161*TAN(A61  •PI()/180)"4-EXP(-MIN(A78.200)) 

isotropic  form  of  pert  theory 

=  10'LOG10(A79j 

pert,  theory  SS 

=  RETURN(MAX{A80.  200)) 

return  perturbation  theory  SS 

=  158"(A63Mreq'{1/3))"(-0.58) 

beta 

=3.3'A82-LOG10(A61/30) 

=42  4‘LOG10(A82) 

=  A83A84t2.6 

Chapman-Hams  formula 

=  RETURN(MAX(A85,-200)) 

return  CH  SS 

Fig.  10  -  F.xcel  spreadsheet  macro  listing  of  surface  scatter  algorithm 
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Harris  formula  does  a  poor  job  of  predicting  surface  scattering  strengths  for  many  combinations 
of  surface  conditions  and  frequencies,  though  there  are  certainly  regimes  where  it  appears  to  do 
a  good  job.  For  those  regimes  where  Chapman-Karris  is  inadequate,  even  a  simple  perturbation 
theory  calculation  such  as  that  contained  in  this  algorithm  does  a  better  job  of  predicting  scat¬ 
tering  strengths.  The  existence  of  a  transition  region  between  air-water  interface  scattering  and 
(presumably)  scattering  from  subsurface  bubbles  is  also  a  reasonable  thing  to  expect;  while  the 
boundaries  we  have  drawn  are  only  approximations,  there  is  considerable  evidence  to  show  that 
they  are  approximately  correct. 

On  the  other  hand,  the  surface  scattering  problem  is  clearly  more  complex  than  this  algorithm 
suggests.  The  principal  objection  to  the  algorithm  is  that  it  is  unlikely  that  surface  scattering 
strengths  can  be  predicted  accurately  on  the  basis  of  a  single  environmental  parameter,  and  that 
even  if  it  could,  that  parameter  would  probably  not  be  instantaneous  wind  speed  measured  at  any 
height.  It  will  be  necessary  in  future  work  to  identify  how  factors  such  as  wind  stress,  wind  history, 
sea  state,  and  water  temperature  affect  surface  scattering  so  that  a  more  complete  model  can  be 
constructed. 
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